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Abstract 

We study a recently proposed mechanism to solve the hierarchy problem in the context of the 
landscape, where the solution of the hierarchy problem is connected to the requirement of having 
baryons in our universe via Electroweak Baryogenesis. The phase transition is triggered by the 
fermion condensation of a new gauge sector which becomes strong at a scale A determined by 
dimensional transmutation, and it is mediated to the standard model by a new singlet field. In a 
"friendly" neighborhood of the landscape, where only the relevant operators are "scanned" among 
the vacua, baryogenesis is effective only if the higgs mass rrih is comparable to this low scale A, 
forcing rrih ~ A, and solving the hierarchy problem. A new CP violating phase is needed coupling 
the new singlet and the higgs field to new matter fields. We study the constraints on this model 
given by baryogenesis and by the electron electric dipole moment (EDM), and we briefly comment 
on gauge coupling unification and on dark matter relic abundance. We find that next generation 
experiments on the EDM will be sensitive to essentially the entire viable region of the parameter 
space, so that absence of a signal would effectively rule out the model. 



1 Introduction 



The most problematic characteristic of the Standard Model (SM) seems to be the smallness of its 
superrenormalizable couplings, the higgs mass m^, and the cosmological constant A,,, with respect to 
the apparent cutoff of the theory, the Planck mass, Mpi. These give rise respectively to the hierarchy, 
and cosmological constant problems. According to the naturalness hypothesis, the smallness of these 
operators should be understood in terms of some dynamical mechanism, and this has been a driving 
motivation in the last two decades in the high energy theory community. 
In the last few years, several things have changed. 

Concerning the hierarchy problem, experimental investigation has shown that already many of 
the theoretically most attractive possibilities for the stability of the weak scale, such as supersym- 
metry, begins to be disfavored, or present at least some fine tuning issues 

Concerning the cosmological constant problem, there was the hope that some symmetry con- 
nected to quantum gravity would have forced the cosmological constant to be zero. However, first, 
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cosmological observations have shown evidence for a non zero cosmological constant in our present 
universe 015; second, from the string theory point of view, two main things have occurred: on the 
one hand, consistent solutions with a non zero cosmological constant have been found 001311111]) 
and, on the other hand, it has been becoming more and more clear that string theory appear to 
have a huge landscape of vacua, each one with different low energy parameters [inilT^ . 

If the landscape is revealed to be real, it would force a big change in the way physics has to be 
done, and some deep questions may find a complete new answer. In particular, it is conceivable that 
some characteristics of our low energy theory are just accidental consequences of the vacua in which 
we happen to be. This is a very big step from the kind of physics we are used to, and its consequences 
have been explored recently in There , it is shown that the presence of a landscape offers a new 
set of tools to address old questions regarding the low energy effective theory of our universe. On one 
hand, there are statistical arguments, according to which we might explore the statistically favored 
characteristics for a vacuum in the landscape [I21IiniIIlllISlIinilIZl[IHlIiniEniEIlE2]- On the 
other hand, there are also some selection rules due to anthropic reasoning, which we might have to 
impose on the vacuum, in order for an observer to be present in it. Pioneering in this, is Weinberg's 
"structure principle" > which predicts the right order of magnitude for the cosmological constant 
starting from the requirement that structures should have had time to form in our universe in order 
for life to be present in it. 

At this point, it is necessary to speak about predictivity for theories formulated with a lot of 
different vacua, as it occurs in the landscape. There are two important points that greatly affect 
the predictivity of the theory. First, we must know how the characteristic of the low energy theory 
change as we scan among the vacua. Second, we must know also how the probability of populate such 
a vacuum changes among the vacua, also considering the influence of the cosmological evolution. 
These are very deep questions whose answer in general requires a full knowledge of the UV theory, 
and we do not address them here. However, there is still something we can do. In fact, as it was 
pointed out in theories with a large number of vacua can be described in an effective field theory 
approach. In such a study, it was shown that it is very natural for parameters not to effectively 
scan in the landscape, unless their average value in the landscape is null. So, it is reasonable to 
assume that some parameters in the landscape do scan, and some do not. Among the parameters 
of the theory, the relevant operators have a particular strong impact on the infrared properties of 
the theory. But exactly because of this, if they are the only ones to scan, their value can be fixed 
by environmental arguments. This is however true only if the marginal couplings do not scan. For 
this reason, a neighborhood of the landscape in which only the relevant parameters are scanning, is 
called a "friendly neighborhood", because it allows to fix the value of the relevant couplings with 
environmental reasoning, while the marginal parameters do not scan, and so they are analyzed 
with the usual instruments of physics, judging them on the basis of their simplicity and of their 
correlations. We do not know if the true string theory landscape have these properties. However, 
there is some phenomenological motivation to expect that this could well be. In fact, in the physics 
connected to the standard model, we have been able to successfully address question concerning the 
marginal parameters with dynamical reasoning, deducing some striking features as gauge coupling 
unification, chiral symmetry breaking, or the weakly interacting dark matter, while, on the contrary, 
we have been having big troubles concerning the problems connected to the relevant coupling in the 
standard model, the cosmological constant, and the higgs scale, both of which arise large fine tuning 
issues. In this paper, we want to address the hierarchy problem of the electroweak scale assuming 
we are living in such a "friendly neighborhood" of the landscape. Because of this, we concentrate a 
little more on the predictivity of a phenomenological model based on this. As it is clear, predictivity 
in the landscape is very much enhanced if there is an infrared fragile feature which must be 
realized in the universe in order for it not to be lethal, and which is however difficult to realize in 
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the scanned landscape. Exactly the fact that a necessary fragile feature is difficult to realize gives a 
lot of constraints on the possible vacua in which we should be, or, in other words, on the parameters 
of the theory, making then the model predictive. An example of a fragile feature is the presence 
of structures in the universe. According to Weinberg's argument |23j . if we require the fact that in 
the universe there should be structures, than the value of the cosmological constant is very much 
tighten to a value close to the observed one. The presence of structures can then also be used to 
anthropically require the lightness of the higgs field in the case dark matter particles receive mass 
from electroweak symmetry breaking |1H I24j . 

In this paper, we concentrate on another fragile feature which we consider necessary for the 
development of any sort of possible life in the universe, and which is a necessary completion on 
Weinberg's structure principle. In the absence of baryons, the dark matter would just form virilized 
structures, and not clumped structures, which instead are necessary for the development of life. 
We can construct a model where baryons are explicitly a fragile feature. Then, we can connect 
the solutions of the problem of baryons to the hierarchy problem with a simple mechanism which 
imposes a low energy electroweak scale in order for baryons to be present in the universe. This is the 
line of thought implemented in which is the model on which we focus in this paper. As we will 
explain in detail in the next section, the mechanism is based on the fact that baryogenesis is possible 
in the model only if the electroweak scale is close to a hierarchical small scale. The hierarchical 
small scale is naturally introduced setting it equal to the scale at which a new gauge sector of the 
theory becomes strong through dimensional transmutation. This mechanism naturally provides a 
small scale. In order to successfully implement baryogenesis, we will also need to add some more 
CP violation. This will force the presence of new particles and new couplings. The model becomes 
also predictive: the lightest of these particle will naturally be a weakly interacting particle at the 
weak scale, and so it will be a very good candidate for dark matter; further, the new CP violating 
terms will lead to the prediction of a strong electron electric dipole. 

The main purpose of this paper is to investigate in detail the mechanism in which baryogenesis 
occurs in this model, and the consequent predictions of the model. In particular, we will find 
that next generation experiments of the electron Electric Dipole Moment (EDM), together with the 
turning on of LHC, are going to explore the entire viable region of the parameter space, constraining 
it to such a peculiar region of the parameter space, so that absence of a signal would result in ruling 
out the model. 

The paper is organized as follows: in sec. II, we explain in detail the model; in sec. Ill we study 
the amount of baryons produced; in sec. IV we determine the electron EDM; in sec. V we briefly 
comment on Dark Matter and Gauge Coupling Unification ; and in sec. VI, we draw our conclusions. 

2 Hierarchy from Baryogenesis 

In this section, following ^j, we show how we can connect the electroweak scale to a scale expo- 
nentially smaller than the cutoff. 

As the presence of baryons is a necessary condition for the formation of clumped structures, it 
is naturally to require that the vacuum in which we are should allow for the formation of a net 
baryon number in the universe. This is not a trivial requirement. In fact, in the early hot phase of 
the universe, before the electroweak phase transition, baryon number violating interactions through 
weak sphalerons were unsuppressed, with a rate given approximately by '^5^: 

r^s = GkaiT (1) 

where A; ~ 20, with the consequence of erasing all precedently generated baryon asymmetry (if 
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no other macroscopic charges are present). However, at the electroweak phase transition, ah the 
Sakharov's necessary conditions [20] for generating a baryon asymmetry are satisfied, and so it is 
possible in principle to generate a net baryon number at the electroweak phase transition, through 
a process known in the literature as electroweak baryogenesis (see [271 128| for two nice reviews). 
However, it is known that the SM electroweak phase transition can not, if unmodified, generate 
the right amount of baryons for two separate reasons. On one hand, CP violating interactions are 
insufficient |29|, and, on the other hand, the phase transition is not enough first order for preventing 
weak sphalerons to be active also after the phase transition [201 • Since this last point is very 
important for our discussion, let us review it in detail. 

Let us suppose in some early phase of the universe we have produced some initial baryon number 
B and lepton number L, and no other macroscopic charge. In particular B — L = 0. The quantum 
number B + L is anomalous, and the equation for the abundance of particles carrying B + L charge 
is given by: 

where ub+l is the number of baryons and leptons per unit volume, s is the entropy density ~ T^, 

Mtt w f (r) 1 v(T) 

and r is related to the sphaleron rate [^Tj, F ~ ^NfTe ^s2'^V2t'^ where is the temperature 
dependent higgs vev {v{Tq) = 246 GeV as measured in our vacuum), Nf is the number of fermionic 
families, and §2 is the SU(2) weak coupling. The reason why reactions destroying baryons are faster 
than reactions creating them is due to the fact that the relative reaction rate goes as: 

^ - e-^/ (3) 

where A/ is the difference in free energy, and T± are the sphaleron rates in the two directions. Now, 
it is easy to see that the free energy grows with the chemical potential fiB, which then grows with 
the number of baryons ub- In the limit of small difference, we then get eq.©. 
This differential equation can be integrated to give: 

rnB+L^ Expf-^4^e-(^^)l (4) 

V S /final V S /initial \ v{Tc) J 

where is the critical density, and where g^: is the number of effective degrees of freedom (~ 55). 

If in the universe there is no macroscopic lepton number, than clearly at present time we would 
have no baryons left, unless: 



9* v{Tc) 

which roughly implies the constraint: 



. 92 V2Tc } <l (5) 



>1 (6) 



This is the so called " baryons wash out" , and the origin of the requirement that the electroweak 
phase transition should be strongly first order. 

Note that this is the same condition we would get if we required sphaleron interactions not to 
be in thermal equilibrium at the end of the phase transition: 
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implies roughly: 

47r v(Tc) 

Tc e «2 V2T^ < ^ (8) 

and so 

^>l (9) 

Also note that the requirement for the sphalerons not to be in thermal equilibrium already just after 
the phase transition is necessary, as otherwise we would have a baryon symmetric universe, which 
leads to a far too small residual relic density of baryons. 

Later on, when we shall study electroweak baryogenesis, we shall get a number for the baryon 
number. In order to get then the baryon number at, let us say. Big Bang Nucleosynthesis (BBN), 
we need to multiply that number by the factor in eq.(@J). In the next sections, we shall consider 
the requirement in eq.® fulfilled, and we will consider completely negligible the wash out from the 
sphalerons coming from eq.ljl]). 

Now, let us go back to the higgs phase transition, and let us assume that in the neighborhood 
of the landscape in which we are, all the high energy mechanisms for producing baryons have been 
shut down. This is easy to imagine if, for example, the reheating temperature is smaller than the 
GUT scale. We are then left with the only mechanism of electroweak baryogenesis. From the former 
discussion, it appears clear that there should be some physics beyond the SM to help to make the 
phase transition strong enough. 

We may achieve this by coupling the higgs field to a singlet S, with the potential equal to: 

V = XS^ + \h{h)h - XS'^ f + ^ml^^ + mlh)h (10) 

where we have assumed a symmetry S — > —S. We can couple this field to two fermions ^ , ^'^^ which 
are charged under a non-Abelian gauge group through the interaction 

ksS^^" (11) 

In order to preserve the symmetry S —S, we give the fields ^, charge i. We can then assume 
that this sector undergos confinement and chiral symmetry breaking at its QCD scale determined 
by dimensional transmutation 

< ^^r^ >r^ A3 (12) 

which is naturally exponentially smaller than the cutoff of the theory. We assume that this phase 
transition is first order, so that departure from thermal equilibrium is guaranteed. 

Now, following our discussion in the introduction, suppose that, scanning in the landscape, the 
only parameters which are effectively scanned are the relevant couplings nih and ms- If then we 
must have baryons in our universe so that clumped structures can form, than we need to be in the 
vacuum in which these two parameters allow for a strong enough first order phase transition in the 
electroweak sector. So, we have to require that this phase transition triggers the electroweak phase 
transition. It is clear that this can only be if it triggers a phase transition in the S field, which 
is possible then only if ms is of the order of A. Finally, the phase transition in S can trigger a 
strong first order phase transition in the higgs field only if again mh ~ 7715 ~ A (for a more detailed 
discussion, see next subsection). So, summarizing, we see that, the requirement of having baryons 
in the universe forces the higgs mass to be exponentially smaller than the cutoff, solving in this way 
the hierarchy problem. 

In order to produce baryons, we still need to improve the CP violating interactions, that in the 
SM are not strong enough. We can minimally extend the introduced model to include a singlet s and 



5 



2 SU(2) doublets ^'±, with hypercharge ±1/2 (notice that they have the same quantum numbers 
as higgsinos in the Minimal Supersymmetric Standard Model (MSSM)), with the following Yukawa 
couplings: 

kSss + k'S'^+'^^+gh^'^+s + g'h'^-S (13) 
There is then a reparametrization invariant CP violating phase: 

e = Bxg{kk'g*g'*) (14) 

The mass terms for this new fields can be prohibited giving proper charges to the fields 
This implies that, since at the electrowcak phase transition these new fermions get a mass of order 
of the electroweak scale, and also because of the fact that the lightest fermion is stable, we actually 
have a nice candidate for dark matter. This last point is a connection between the higgs mass, 
which, up to this point, we have just assured to be exponentially smaller than the cutoff, and the 
weak scale, but this connection will come out quite naturally later. So, if this model happens to 
describe our universe, what we should see at LHC should be the higgs, the two new singlets S and 
s, and the two new doublets 

Since the model is particularly minimal, it is interesting to explore the possibility for generating 
the baryon number of the universe in more detail. Before doing so, however, let us see in more detail 
how the vevs of the higgs and of the singlet S are changed by the phase transition. 



2.1 Phase Transition more in detail 

Before going on, here we show more in detail how the requirement of having a strong first order 
phase transition leads to have ms ~ ruh ~ A. 

In unitary gauge, the equations to minimize the potential are (from here on, we mean by S also 
the vev of the field S\ the meaning will be clear by the contest): 

- 4.\h\S{— - XS^) + mis + ks^'Q{Tc -T) = (15) 

2Xn{^-~XS')^ + 2ml^ = (16) 

where Tc is the critical temperature. The first equation comes from the derivative with respect to 
S, and we will refer to it as S equation, while for the other we will use the name h equation. 

We first consider the case m| > and m| > 0. Before the phase transition, we have the 
minimum at the symmetric vacua: 

5 = 0, v = 0, for T > (17) 

For T < Tc, the minimum conditions change, and we can not solve them analytically. We can 
nevertheless draw some important conclusions. Let us consider the minimum equation for S, which 
is the only equation which changes. Let us first consider the case ms S> A. In this case, we can 
consistently neglect the cubic terms in the S equation, to get: 

Sc.-^^ (18) 

{ml-2XXhV^) 



Then, if m| » 2XXhv'^, we have: 



S = (19) 



mg 
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v = (20) 
while the other solutions are still unphysical, as: 



9 2 

V m 



(fc|A6) 



j_^_:^ + aa::;^ <o (21) 

If m| <C 2AA/jU^, then: 



2 ' '■ m| 



5:.--^ (22) 
2AAfti;2 ^ ^ 



For the higgs, the non null solution is: 



which has some relevant effect for the electroweak phase transition only if nih ~ A. But in that case 
-y ~ A, and from m| ^ 2AA/jV^ we get that m| ^ A^, in contradiction with our initial assumption. 

Then, if "^s ^^^''^ ^ 1, in the S equation, we can consistently consider just the cubic term, to get: 



Sc^ ^ (24) 

4i/3(A + AA,,)V3 ^ ' 

then, the non null solution for v'^ becomes: 

— ~ — + A — M- = (25) 

2 \h \AV^{X + XXn){X + XXh)) 

which has some relevant effect only if m\ -C A^. However, in this case the condition ^^"^^^ ^ 1 
would imply ^ <C 1, again in contradiction with our assumptions. 

So, in order to have some effect on the higgs phase transition, we are left with the only possibility 
of having ms <C A. 

Restricting to this, consistently, we can neglect the linear term in the S equation, to have: 

4A5^ - 4A,,A5 - AS^^ = -kgk^ (26) 

2A^ (y - A^^^ = -2ml (27) 
which implies the following equation for S: 

AXS^ + 4Am^5 = -ksK^ (28) 

For ruh S> A, we have: 

and the non null solution for v is still not physical: 



So, finally, we have that in order to have a strong first order phase transition triggered by the new 
sector, we are forced to have ruh <C A, which is what we wanted to show. In detail, for m/j <C A, we 
implement the following phase transition 



Now, the final step is to show that, not only nih <C A and ms <C A, but that actually nih ~ A 
and ms ~ A. The argument for this is that, scanning in the landscape, it will be generically much 
more difficult to encounter light scalar masses, as they are fine tuned. So, in the anthropically 
allowed range, our parameters shall most probably be in the upper part of the allowed range. So, 
we conclude that this model predicts nih ms ~ A, as we wanted to show. 

This phase transition must satisfy the requirement that the sphalerons are ineffective (^^^^ ^ 1); 
after the phase transition, which occurs at, roughly, the critical temperature Tc ~ A. So: 

^.qi|)!^!>, (3.) 

which can clearly be satisfied for some choices of the couplings. In order to better understand the 
natural values of the ratio v{Tc)/Tc in terms of the scalar couplings, it is worth to notice that the 

coupling A appears in the lagrangian as always multiplied by \h- So, the coupling which naturally 
tends to be equal to the other ones in the lagrangian is Ag = A/jA. With this redefinition, we get the 
constraint: 

n UJ V A J 

which can clearly be satisfied with some choice of the parameters. It is also worth to write the ratio 
of the vevs of S and h after the phase transition: 

We can also analyze the case in which m| < 0. In this case, the electroweak phase transition 
would occur before the actual strong sector phase transition. However, since we know that in the 
SM the phase transition is neither enough first order, nor enough CP violating, we still need, in 
order to have baryon formation, to require that the strong sector phase transition triggers a phase 
transition in the higgs sector. It is then easy to see that all the former discussion still applies with 
tiny changes, and wc get the same condition m|. ~ A^. There is a further check to make, 

though, which is due to the fact that, in order for baryogenesis to occur, wc need an unsuppressed 
sphaleron rate in the exterior of the bubble. This translates in the requirement, for < 0: 

1^ « 1 (35) 

which can be satisfied in some vacua. 

In the next of this paper, we will always assume that these conditions are satisfied, and the 
sphalerons are suppressed in the broken phase. 

Now, we are ready to treat baryogenesis in detail. 



> 1 (33) 
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3 Electroweak Baryogenesis 



During the electroweak phase transition, we have ah the necessary conditions to fulfill baryogenesis 
|27l 128] . We have departure from thermal equilibrium because of the phase transition; we have CP 
and C violation because of the CP and C violating interactions; and finally we have baryon violation 
because of the unsuppressed sphaleron rate. The sphaleron rate per unit time per unit volume is in 
fact unsuppressed in the unbroken phase (see eq.Q). 

There are various different effects that contribute to the final production of baryons. For example, 
CP violation can be due just to some CP violating Yukawa coupling in the mass matrix, or it can 
be mainly due to the fact that, in the presence of the wall, the mass matrix is diagonalized with 
space-time dependent rotation matrixes, which induce CP violation. Further, CP violation can be 
accounted for by some time dependent effective coupling in some interaction terms. Baryon number 
production as well can be treated differently. In the contest of electroweak baryogenesis which we 
are dealing with, there are mainly two ways of approaching the problem, one in which the baryon 
production occurs locally where the CP violation is taking place, so called local baryogenesis, or one 
in which it occurs well in the exterior of the bubble, in the unbroken phase, in the so called non 
local baryogenesis. At the current status of the art, it appears that the non local baryogenesis is 
the dominant effect, at least for not very large velocity of the wall v^, which however is believed to 
be not large because the interactions with the plasma tend to slow down the wall considerably j31j . 
and we concentrate on this case (see Appendix A for a brief treatment of baryogenesis in the fast 
wall approximation). 

In order to compute the produced baryon abundance, we follow a semiclassical method developed 
in |S2]- A method based on the quantum Boltzmann equation and the closed time path integral 
formalism was developed in [33], making a part of the method more precise. However, the corrections 
given by applying this procedure to our case can be expected to be in general not very important 
once compared to the uncertainties associated to our poor knowledge of certain parameters of the 
electroweak phase transition, as it will become clear later, which make the computed final baryon 
abundance reliable only to approximately one order of magnitude |32j. Further, it is nice to note 
that our general conclusions will be quite robust under our estimated uncertainty in the computation 
of the baryon abundance. For these reasons, the method in ^] represents for us the right mixture 
between accuracy and simplicity which is in the scope of our paper. 

Since the method is quite contorted, let us see immediately where the basic ingredients for 
baryogenesis are. Departure from thermal equilibrium obviously occurs because of the crossing of 
the wall. C violation occurs because of the V — A nature of the interactions. CP violation occurs 
because the CP violating phases in the mass matrix are rotated away at two different points by 
two different unitary matrix such that U{x)^U{x + dx) 1. There could be other sources of CP 
violation, but, in our case, this is the dominant one. Baryon production occurs instead well in the 
exterior of the wall, in the so called non-local baryogenesis approach, where the weak sphaleron rate 
is unsuppressed. 

Let us anticipate the general logic. The calculation is naturally split into two part. In the first 
part, we compute the sources for CP violating charges which are due to the CP violating interactions 
of the particles with the incoming wall. This calculation will be done restricting ourself to the vicinity 
of the wall, and solving a set of coupled Dirac-Majorana equations to determine the transmission 
and reflection coefficients of the particles in the thermal bath when they hit the wall, which are 
different for particles and antiparticles. In the wall rest frame, the wall is perceived as a space-time 
dependent mass term. This will give rise to a CP violating current. The non null divergence of this 
current will be the input of the second part of the calculation. In this second part, we shall move 
to a larger scale, and describe the plasma in a fluid approximation, where we shall study effective 
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diffusion equations. Tlie key observation is tliat, once the charges have diffused in the unbroken 
phase, thermal equilibrium of the sphalerons will force a net baryon number. In fact, in the presence 
of SU(2) not neutral charges, the equilibrium value of the baryon number is not zero: 

(B + L)e, = c,Qi (36) 

i 

where Ci are coefficients which depends on the different charges . Once produced, the baryon number 
will diffuse back in the broken phase, where, due to the suppression of the sphaleron rate, it will be 
practically conserved up to the present epoch. This will end our calculation. 

In the next two subsections, we proceed to the two parts of the outlined calculation. 

3.1 CP violation sources 

We begin by computing the source for the CP violating charges, following [^. We restrict to the 
region very close the the wall, so that the wall can be considered flat, and we can approximately 
consider the problem as one dimensional. We consider a set of particles with mass matrix M(z) 
where z is the coordinate transverse to the wall, moving , in the rest frame of the wall, with energy- 
momentum E, k. Taken zq as the last scattering point, these particles will propagate freely for a 
mean free time Tj, when they will rescatter at the point ZQ + TiV, where v is the velocity perpendicular 
to the wall, ktr/E. Now, because of the space dependent CP violating mass matrix, these particles 
will effectively scatter, and the probability of being transmitted and reflected will be different for 
particles and antiparticles. This will create a current for some charges, whose divergence will then 
be the source term in the diffusion equations we shall deal with in the next section. The effect will 
be particularly large for charges which are explicitly violated by the presence of the mass matrix, 
and we shall restrict to them. 

We introduce J± as the average current resulting from particles moving towards the positive 
and the negative z direction, between zq and zq + A, where A = tv, and r is the coherence time of 
the particles due to interactions with the plasma jHSj- J± are the CP violating currents associated 
with each layer of thickness A. J_|_ receives, for example, contributions from particles originating 
from the thermal bath at zq with velocity v, and propagating until zo + A, as well as from particles 
originating at zq + A with velocity —v, and being reflected back at zq + A. The formula for J± is 
given by: 

j+ = (tt (t^qt - Ttgr)) - Tt (p,,+a [r^QR - A^QR))) (l,0,0,t5) (37) 

J_ = (Tr (^R^QR - R^QRJ ) - Tr (p,„+A (f^QT - T^Qt) ) ) (1, 0, 0, -r;) (38) 

where R{R) and T[T) are reflection and transmission matrices of particles (antiparticles) produced at 
Zq with probability pz^, evolving towards positive z; while T and R are the correspondent quantities 
for particles produced at + A with probability Pz^+A and evolving towards negative z; v is the 
group velocity perpendicular to the wall at the point zq, and v is the same quantity at zq + A; Q is 
the operator correspondent to the chosen charge; and the trace is taken over all the relevant degrees 
of freedom and averaged over location zq within a layer of thickness A. When boosted in the plasma 
rest frame, these currents will become the building block to construct the CP violating source for the 
charges that, diffusing into the unbroken phase, will let the production of baryon number possible. 

Now, consider a small volume of the plasma, in the plasma rest frame. As the wall crosses 
it, it leaves a current density equal to (7+ + J-)niasma ^very time interval r, where the subscript 
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plasma refers to the quantity boosted in the plasma frame. So, at a time t, the total current density 
accumulated will be given by: 



dt'kj^{x,t') + J^{x,t'))Zasrna (39) 



T 



where tr is the relaxation time due to plasma interaction. Prom this, the rate of change of the 
charge Q per unit time is given by: 

7Q(f,t)= d^s^= (40) 
1 1 

-(J+(f,t) + J__(f,t))J;„- -(J+(x,t-r/j) + J_(x,t-rR))°;„ 

— / dzdz{J+ + J-Yplasma 

Since in the SM CP violation is very small, and already proved to be not enough to account for 
the baryon number of the universe, we can clearly concentrate in the sector of the neutral particles 
of the new theory, where the mass matrix is given by: 



M{z) 



The transmission and reflection coefficients can be found by solving the free coupled Dirac-Majorana 
equations for these particles with the mass matrix given in M{z). We can solve this by a method 
developed in [32], in a perturbative expansion in mass insertion. As it is explained in |321 I34j . 
the small parameter in the expansion is mA, where m is the typical mass of the particles, in the 
case wA < 1, or m/w in the case wA > 1. In both cases, the expansion parameter is smaller 
than one. This can be understood noticing the analogy of our system with the scattering off a 
diffracting medium with a step potential of order m. In that case, reflection and transmission are 
comparable (and this is the only case in which we produce a net CP violating charge) only if the 
wave packet penetrates coherently over a distance of order 1/m, and has few oscillation over that 
distance. Suppression of the reflection occurs both if niA <^ 1 and ii ni/w <^ 1. In the first case, 
this is because only a layer of thickness A contributes to the coherent reconstruction of the reflected 
wave, while in the second case, because fast oscillations tend to attenuate the reconstruction of the 
reflected wave. Up to sixth order in the mass insertion, we get: 

T= 1- [ dzi /" 'dz2M2Mi*e2™(^i-^2) + (41) 
Jo Jo 

+ [ dzi r dz2 f (iz3 / 'dz4M4M3*M2 Ml* e^^^^i -^2+^3-^4) + 

Jo Jo Jz2 Jo 

- I dzi ( ' dz2 I dz3 [ ' dzi I dz5 I 'dz6M6M5*M4M3*M2Mi*e2™(^i-^2+23-^4+25-26) + 

Jo Jo Jz2 Jo J Z4, Jo 



+ ... 
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f= I- dzi dz2M^Mie-^''"^'^-'^^ + 

Jo Jzi 

+ [ dzi ( dz2 [ ' dz3 [ dz4M4*M3M2*Mie-2*™(^^-^2+^3-^4) + 

Jo Jzi Jo Jz3 



pA pA pz2 rA pz4 pA 

/ dzi / dz2 / dz3 / dz4 / dz^ / dz6M^M5MlM3M^Mie-'^''"^''^"''^+''^~'^+'""'''^ 

Jo Jzi Jo Jz3 Jo Jz5 



R= -I dziM*e2^"'(^i) + (43) 







+ / dzi [ ' dz2 [ dz3M3*M2Mi*e2^"'(^i-^2+^3) + 

Jo Jo Jz2 

- I dzi [ ' dz2 [ dz^ r dzi [ (iz5M*M4M3*M2Mi*e2™(^i-^2+^3-24+^5) + 

Jo Jo Jz2 Jo J Zi 



+ ... 



[■A 

R= / (iziMie-2™(^i) + (44) 
Jo 

- [ dzi [ dz2 [ ' (iz3M3M2*Mie-2™(^i-^2+^3) + 

rA r-A (■Z2 l-A r-Zi 

+ dzi dz2 / dz3 / dz4 / (i25M5M4*M3M2*Mie-2™(^i-^2+23-^4+^5) ^ 
JO ^0 ^^3 Jo 



where Mi = M(zj).The analogous quantities for the antiparticles are obtained replacing M M* 
in all the former formulas. 

We also need to have the density matrices and p^o+A-We can choose these densities as 
describing thermal equilibrium densities in eigenstates of the unbroken phase: 

Pzo =I)iag {ns{E,v),n^_^{E,v),n^^{E,v)) (45) 

where n.(E, v) is the Fermi-dirac distribution, boosted in the wall frame: 

n = ,^ , (46) 

e +1 

and Pzo+A = Pzo{v The motivation of this is that particles are produced in interaction 

eigenstates which differ from mass eigenstates by a unitary rotation; ignoring this, amounts at 
ignoring small corrections of order {M{z)/Tc)'^ . The choice of thermal distribution is particularly 
good in the small velocity regime, in which we have restricted, where the non thermal contribution 
is of order f^, and it induces corrections of order v'^ in the final baryon density |32j . 

Finally, we have to consider the charges which can play a role in generating the baryon number. 
When choosing such charges, one has to consider that the most important charges are those which 
are approximately conserved in the unbroken phase, as these are the ones which can efficiently diffuse 
in the unbroken phase, and induce a large generation of baryon number. Keeping this in mind, it 
is easy to see that the only relevant charge in our model is the "higgs number" charge, which in 
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the same basis in which we expressed the mass matrix, for the new CP violating sector particles, is 
given by: 

Q,, = Diag (0,0, 1,-1) (47) 

The name "higgs number" just comes from the fact that the fields ^± have the same quantum 
numbers as higgsinos in the MSSM. Now, we can substitute in the formulas for J+ and J_. In 
order to keep some analytical expression, we decide to do a derivative expansion M{z) = M(zo) + 
M'{zq){z — zq) + 0{t/w)'^ and v = v + 0{t /w)"^. This expansion is justified in the parameter range 
T <^ w, which is expected to be approximately fulfilled 

Unlike in the case of the analogous calculation for the MSSM |32] where the leading effect occurs 
at fourth order in mass insertion, here the leading contribution occurs at sixth order. We will explain 
later the reason of this. For the moment, we get: 



J+= (1,0,0,^;) X (48) 

f-A f-zi f-A j-z-i f-A f-zs 

+4 / dzi / dz2 / dzs / dz4 / dz^ / dzQsm{2w{zi - Z2 + z^ - Z4^ + z^ - zq)) 

Jo Jo Jz2 Jo Jz4 Jo 

Im {p,,QhMeM^M^MlM2Ml) + 

pA pzi pA pz'^ pA pz^y 

—4 / dzi I dz2 / dz3 / dz4 / dz^ / dzg sm{2w{zi — Z2 — z^ + Z4 — z^ + zq)) 
Jo Jo Jo Jo JzA Jo 



Im {p^^QhU^M^M^M^MiMl) + 
(■A pA pA pz3 



t-lX t-lX pIX i-z-i i-IX pzs 

+4 / dzi / dz2 / dz3 / dz4 / dz^ / dze sm{2w{zi - Z2 + Z3 - z^ + zrj - zq)) 
Jo Jo Jzo Jo JzA Jo 



'22 JO J Z4, 

Im {p,,+AQhMlM6M;MiM^M2) 



J_ = (1,0,0,-v) X (49) 

pA pA pz2 pA pz4 pA 

—4 / dzi / dz2 / dz^ / dz4 / dz^ / dzQ sin(2w(2:i — Z2 + z-i — Z4 + z^ — zq)) 

Jo Jzi Jo Jz:i Jo Jz5 

Im {p,,+AQhM6M*M4M*M2M*) + 

i-A j-A j-A i-A rzi pA 

+4 / dzi / dz2 / dz^ / dz4 / dz^ I dze sm{2w{zi - Z2 - Z3 + Z4 - Z5 + zq)) 

Jo Jzi Jo Jz3 Jo Jz5 

lm{p,,+AQhMiM;MeM;M4Ml) + 

pA pA pz2 rA pz4 pA 

—4 / dzi / dz2 / dz^ / dz4 / dz^ / dzQ sin(2tt;(2:i — Z2 + z-i — Z4 + z^ — zq)) 

Jo Jo Jo Jz3 Jo Jzs 

Im {p,,QhMlM(,MlM4M*^M2) 

where the z dependence in each mass matrix Mi is to be understood at linearized level. 

We can substitute the results in eq.(|ir?|). to get an expression for the higgs charge source. In 
reality, if we take the relaxation time large enough, and if we keep performing a derivative expansion, 
only the first term in eq. (|4()j) is relevant. At first order in the wall velocity, we get: 

7Q = J X (50) 



i=ip,s 



d^kfiwiA)^ e-^'/^- Ei 



27r3 w'^ \i + e.E^IT.)2T^ 
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where ^7 is a CP violating invariant: 



J=\gg'kk'{g^-g'^)sm{e) (51) 
and fi{wA) is: 

hi^^) = ''^'fjr'"\ {S'{zo){9 - 2Wa2 - Uw^A^Mzo) - 9v'{zo)S{zo) (52) 
-2w^A(llt>'(zo)A(3 + w^A'^) + 

3(5 + 2w'^A^)v{zo))S{zo) + 6cos{6wA){-S'{zo)v{zo) + v'{zo)Sizo)) + 
6cos(4u'A)(25'(zo)(-l + 3w'^A'^)vizo) + 2v'izo)S{zo) + w"^ A{-3v' (zq) A + t;(zo))5(zo)) + 
3cos(2u;A)(3S"(zo)(l + Gw"^ A'^)v{zo) - 3v'{zo)S{zo) + 2^)"^ A{3v' (zq) A + 4t>(zo))5'(zo)) 
+2w{2S'{zo)A{-3 + Uw^A'^)v{zo) + (5t;'(zo)A(3 + Tw^A^) + 
3(1 + 7w^A^)v{zo))S{zo)) sm{2wA)3w{16S' izo)Av{zo) 

+(-13^;'(zo)A + vizo))Sizo)) sin(4u>A) + 6wA){-S' izo)v{zo) + v'izo)S{zo)) sin(6u'A)))) 

/,(^A) = !!i£^£2l(S'(zo)(15 + Uw'A' - 28w'A^)v{zo) (53) 
-l5v'{zo)S{zo) - 4w^A{v'{zo)A{A8 + iWA^) + 

3(5 + 2w'^A'^)v{zo))S{zo) + l2cos{6wA){S' {zo)v{zo) - v'{zo)S{zo)) + 3cos{AwA){S' (zq) 
(1 - 12w^A'^)v{zo) - v'izo)Sizo) + A{6v' (zq) A + v{zo))S{zoO) + 
6cosi2wA{S'izo){-5 + Uw'^A'^)v{zo) + 5v'{zo)Sizo) +2w^A{5 v'{zo)A + 
Av(zo))S{zo)) + l2w{S'{zo)A{-7 + 6u;2a2)7;(zo) + {v{zo) + 
A(5f'(zo)(2 + 3u'2a2) + 7w'^Av{zo))S{zo)) sin(2?i;A) 

-6w{-9S'izo)Av{zo) + (12?;'(zo)A + v{zo))S{zo)) sin(4u'A) + l2wA{S' {zo)v{zo) 
-v'{zo)S{zo)) s'm{6wA))) 

where zq here is a typical point in the middle of the wall. We can estimate the coherence time as 
r ~ {gwTc)~^, and so we can take, roughly, the typical interval which is valid also in the MSSM '32j: 
15 < tTc < 35, and in the approximation of making the particles massless, which is still consistent 
with our approximations, we can integrate numerically the integral, and fit it linearly in r (which is 
a good approximation in the range of interest), to obtain: 

7Q(f , t) ~ 24t;^7^ {gg'kk'ig^ - g'^) sin(0)) x (54) 
(1 + |(tT,-25)) v^S 



TTr. T} 



{S'v - v'S) 



Even though the reparametrization invariant CP violating phase 9 requires only 4 Yukawa couplings 
to be present in the theory, it is somewhat puzzling that the leading effect appears at sixth order in 
the couplings. The reason is as follows. In the case g = g' , there is an approximate Z2 symmetry 
which exchanges the ^± fields. This is only an approximate symmetry, because ^± differ by their 
hypercharge. However, the hypercharge is considered a subleading effect here, and never comes into 
play in our computation, and so we do not see the breaking of this symmetry. The Z2 symmetry 
has, in the basis we have been using up to now, the matrix form: 
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We clearly see that: 

\M{z), Z2] = 0, [p{z), Z2] = 0, {Q, Z2} = 



(55) 



which then implies that: 

Tr (^p,,TQT^^ = (56) 

and similar for the others terms in J_|_ and J_. So, in the limit in which g = g' , the baryon asymmetry 
should vanish. Then, since each particle needs to have an even number of mass insertions to be 
transmitted or reflected, and since the CP violation requires at least 4 mass insertion, we finally 
obtain that the coupling dependence in the CP violating source must be of the form present in 
ea.(|54|l. From this discussion, it is clear that in the case g = g' the baryon production will be 
heavily suppressed, and, from what we will see later, it will be clear that the couplings in this case 
will have to be so large to necessarily hit a Landau pole at very low energies, making the model 
badly defined. We shall neglect this degeneracy of the couplings for the next of the paper. 
Now, we are ready to begin the second part of the computation. 

3.2 Diffusion Equations 

Here, we begin the second part of the calculation, still following We turn to analyze the system 
at a larger scale, and approximate it to a fluid. We then study the evolution of the CP violating 
charges due to the presence of the sources and of the diffusion effects in the plasma. To this purpose, 
we shall write a set of coupled differential equations which include the effects of diffusion, particle 
number changing reactions, and CP violating source terms, and we shall solve them to find the 
various densities. We shall be interested in the evolution of particles which carry some charges 
which are approximately conserved in the unbroken phase. Near thermal equilibrium, which is a 
good approximation for small velocities, we can approximate the number density as: 

Ui = kiPiT^/6 (57) 

where /Xj is the chemical potential, and ki is a statistical factor which is equal to 2 for each bosonic 
degree of freedom, and 1 for each fermionic. 

The system of differential equations simplifies a lot if we neglect all couplings except the gauge 
couplings, the top quark Yukawa coupling, and the Yukawa couplings in the new CP violating sector. 
From the beginning, we take the interactions mediated from these last ones to be fast with respect 
to the typical timescale of the fluid. We include the effect of strong sphaleron, but neglect the one of 
the weak sphalerons until almost the end of the computation. This allows us to forget about leptons. 
We need only to keep track of the following populations: the top left doublet: Q = {ti + bL)-, the 
right top: T = tR, the Higgs particle plus our new fields ^±:H = /iq + + Strong sphalerons 
will be basically the only process to generate the right bottom quarks B = b^, and the quarks of the 
first two generations Q(i.2)l, Ur, Cr, Sr, Dr. This implies that all these abundance can be expressed 
in terms of the one of B: 

QiL = Q2L = -2Ur = -2Dr = -2Sr = -2Cr = -2B = 2{Q + T) (58) 

The rate of top Yukawa interaction, Higgs violating process, and axial top number violation are 
indicated as ry,r/j,rffj, respectively. We take all the quarks to have the same diffusion equation, 
and the same for the higgs and the ^'s. The charge abundances are then described by the following 



15 



set of differential equations: 

Q= DqV^Q-Ty{Q/kQ-H/kH-T/kT)-Tra{Q/kQ-T/kT) (59) 

-QVss{2Q/kQ - T/kT + 9(Q + T)/kB) 
t= DgV^T-Vyi-Q/kQ + H/kH + T/kT) 

-Tmi-Q/kQ + T/kr) + 3Tss{2Q/kQ - T/kr + HQ + T)/kB) 
H = DhV^H - Tyi-Q/kg + T/kr + H/kn) - ThH/kn + 7q 

We can restrict ourselves to the vicinity of the wall, so that we can neglect the curvature of the 
surface, and assume we can express everything in a variable z = \r + Vwt\. The resulting equations 
of motions become: 

v^Q' = DqQ" -Ty{Q/kQ-H /kn- T/kr)- VmiQ/kq- T/kT) (60) 

-QTss{2Q/kQ - T/kT + 9(Q + T)/kB) 
v^T' = DgT" - Tyi-Q/kg + H/kH + T/kT) 

-Tmi-Q/kQ + T/kT) + 3Tss{2Q/kQ - T/kT + HQ + T)/kB) 
v^H' = DhH" - Ty{-Q/kQ + T/kT + H/ku) - ThH/kf, + Jq 

We now assume that Ty and T^s are very fast, and we develop the result at 0{1/Ty,l/Tss). This 
allows to algebraically express Q and T in terms of H, to get the following relationships: 

^_rr f kgjgkT - kB) \ 

^~ \kH{kB + 9kQ + 9kT)J ^ ' 



r^^_JJ ( kT{2kB+9kQ) 

\kH{kB+9kQ + 9kT), 
and, substituting back, we find the following effective differential equation for H: 



(62) 



v^H' = DH" -TH + j (63) 
where the effective couplings are given by: 

^ ^ Dq{9kQkT - 2kQkB - 2kBkT) + Dhkuj^kQ + 9A:t + fcg) 
9kQkT-2kQkB-2kBkT + kH{9kQ + 9kT + kB) 

^ = ( kn{9kQ+9kT + kB) \ 

' "^\9kQkT-2kQkB-2kBkT + kH{9kQ+9kT + kB) ^ ' 

f^.p ^ .( {9kQ + 9kT + kB) \ 

\T-2kQkB-2kBkT^kH{9kQ+9kT^kB)) ^ ' 



We can estimate the relaxation rates for the higgs number and the axial quark number as |32j : 

(r.„ + r,)~i^|t|i),. 

where At is the top Yukawa coupling, and is the W boson mass, and, in order to keep analytical 
control, we approximate the source term and the relaxation term as step functions: 

7 = 7, w > z>Q (68) 
7 = 0, otherwise 



16 



aoad 

f = f, z>0 (69) 
f = 0, z <0 

For the source term 7, we can take the avaraged value of expression ()54() . However, due to our 
lack of knowledge of the details of the profiles of the fields during the phase transition, we can just 
approximate that expression with: 

7Q = 2Av^^^ {gg'kk'ig^ - g'^) sm{0)) x (70) 
(l + |(rr,-25)) v^S^ 

tTc r3 w 

where we have taken S' ~ S/W and v' ~ v/W, with W the wall width, and we have assumed, as 
expected, that no cancellation is occurring. Here S and v are taken to be of the order they are today. 
In the approximation that D is constant, and with the boundary conditions given by H{zizoo) = 0, 
we have an analytical solution in the unbroken phase j32j : 

H = Ae^""-^!^ (71) 
where, in the limit that Df ^ u^, which is in general applicable, 

2 /"i _ p-'iW^JTiL 



^-^(^l-e-^^ViV^j (72) 

Note that diffusion of the higgs field, and so of the other charges, in the unbroken phase, occurs for 
a distance of order z ~ D /vw 

We now turn on the weak sphaleron rate, which is the responsible for the baryon generation. 
The baryon density follows the following equation of motion: 

VyjpB = Dqp"^ - e{-z)?,Tu,snL{z) (73) 

where we have assumed that the weak sphaleron operates only in the unbroken phase, and where 
is the total number density of left fermions. The solution to this equation is given by, at first 
order in v^: 

PB = -^-^ r d~^^L{-z) (74) 



J —00 

Now, in the approximation in which all the particles in our theory are light, we have: 

kq = 6, kx = 3, ks = 3, kn = 8 (75) 

and in the limit as 00, the resulting baryon abundance is zero P^. This means that 

we have to go to the next order in the strong sphaleron rate expansion. Note also that, with this 
particle content, using the SM quark and higgs diffusion equation Dg ~ G/T^, ~ llO/Tc we 
have D ~ 96/Tc and DTy/v'^ S> so that the assumption that the Yukawa interaction is fast is 

self-consistent, since Ty ~ (27/2)X'^asTc 32 . Finally, we take: 

T^s = QkaljTc, Tss = e/c'-afTc (76) 

where k' is an order one parameter and k 20 25 . To go to next order in the expansion in large 
Tss, we write: 

r._rr( kQ{9kT-kB) \ 

^-''[kn{ks + 9kQ + 9kr))^^^ ^''^ 
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^ - Wiks + 9kQ + 9kr)) + 

Substituting in (I^U)) . we get: 

- 1 — f:^ — ) W{9kQ+9h+k,r) + ^/^^) (^9) 

Using = g + QiL + Qai = 5Q + 4T = (^^^^^) <5q, we get: 

3 f DgH"-v^H' \ 

"^ = ^'^« = "mV J ^ ^ 

Substituting the solution for the Higgs field, we finally get: 

PB _ ( 9ATy,s \ ^\ (8^) 



s \112 s TssJ V D 

where s = {2ir'^ gj 45)Tj^ ~ SST^. Substituting our parameters, we get: 



...... „.(^^)_(l±i^ 



{gg'ig'-g'')kk'sm{e)) 



V \ I b 



It is worth to make a couple of small comments on the parametric dependence of this expression. 
The dependence of the coupling terms, and therefore on the vevs of h and S, was explained in the 
former section. The wall width W has simplified away because the exponential in ea. (|72j) is small 
and can be Taylor expanded. The presence of Tgs in the denominator is due to the particular particle 
content of this model, for which the leading term in the expansion in l/Tgs is zero. The factor 10~^ 
is mainly due to the factor of ~ 10~^ from the entropy density, and the ratio 20 a^/af , times some 
other factors coming from the diffusion terms. For typical values of the wall velocity, we can take 
approximately the SM range: Vy^ ~ 0.05 — 0.3, while the mean free time is r 20/Tc — 30/Tc. With 
these values, the produced baryon number ranges in the regime: 

^ c (2 X 10-^ - 3 X 10-6) [gg'ig' - g")kk' sm{e)) (^)' (|)' (83) 

This is the number which has to be equal to the baryon density at BBN: 

(PR) =9x10-11 (84) 

V S / BBN 

For the moment, let us try do draw some preliminary conclusions on what this does imply on 
the parameters of the model. Clearly, there are still too many parameters which could be varied, so, 
as a beginning, we can set all the Yukawa couplings in the new CP violating sector to be roughly 
equal to each other k ^ k' ^ g ^ g' , and v/Tc ~ 1, we then get: 



(2x 10-^-4x 10-6) (/sin(0)) (^-J -10-^° (85) 
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We then get the following constraint: 

/ (^^J sin(0) ~ 10"* (86) 
It is natural, in this theory, to take the CP violating phase to be of order 1. In this case: 

g g' ^ k k' 0.21 (^^y^' (87) 

From this we see that the assumption of considering the interactions mediated by these Yukawa 
couplings to be fast is justified for a large fraction of the parameter space. These values become 
lower bounds for the couplings if we allow for the CP violating phase to be smaller than order 1 (see 
eq.®). 



4 Electric Dipole Moment 

The same CP violating phase which is responsible for baryogenesis, induces an electric dipole moment 
through the 2-loop diagram shown in figO 




f' f 

Figure 1: 2 loops diagram contributing to fermion EDM, where Xo ^.re the neutral mass eigenstates, 
and L0-^. is the charged one. 

The situation here is much different than in the MSSM, where a CP violating phase in general 
introduces EDM at one loop level, the constraints on which generically force very small CP violating 
phases in the MSSM. It is instead much more similar to the case of split supersymmetry |44l 1451 117] . 
where all the one loop diagrams contributing to the EDM are decoupled. Here the leading diagram 
is at two loops level, and, as we will soon see, it will induce electron EDM naturally just a little 
beyond the present constraints, and on the edge of detection by future experiments. 

The induced EDM is (see [SHI): 



df 



a(r°,r±,r;,) 



a^mf 



2„4 



M, 



E 



dz 



7 Jo 

— / dyy 

7 7o 



Im (OfOf*)g (r°,r±,r/,) , 



dy 
{R- 



yz{y + z/2) 
[z + Rf{z + Ki) 
3Ki)R + 2{Ki + R)y 



AR{Ki - Rf 



+ 



Ki{Ki 



2{Ki 



-2y), Ki 



(90) 
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with 



fl = , + /.. = ^ + -,.*.-^..«.^,.,.^. (91) 

Of = JVJ, of = -Af4,:C« (92) 

where = e~'^ and N'^Mj^^N = diag{m^j,m^2!™x3} with real and positive diagonal elements. 
The plus(minus) sign on the right-hand side of ea. (|88|) corresponds to the fermion / with weak isospin 
+(— )l/2. and /' is its electroweak partner. We mean by to the charged mass eigenstate, and with 
Xi, i = 1,2, 3, the neutral mass eigenstates in order of increasing mass. 

Diagonalizing the 3x3 mass matrix numerically, we see that generically the model predicts EDM 
very close to detection. The induced EDM for some generic parameters are shown in figEl 



de/(e cm) 




100 150 200 300 500 700 1000 1500 2000 



Figure 2: The predicted EDM in this model. We plot the induced electron EDM as a function of 
k'S. The solid line represents the induced EDM with g = g' = 1/2 and kS = 100 GeV, while the 
dashed line represents g = g' = 1/W and kS = 100 GeV, and we take maximal CP violating phase. 
The horizontal line represents the present electron EDM constraint de < 1-7 x 10~^'^e cm at 95% CL 

It is clear then that improvements of the determination of the EDM are going to explore the most 
interesting region of the parameter space. Ongoing and next generation experiments plan to improve 
the EDM sensitivity by several orders of magnitude within a few years. For example, DeMille and 
his Yale group |^ will use the molecule PbO to improve the sensitivity of the electron EDM to 
10~^^ e cm within three years, and possibly to 10"'^^ e cm within five years. Lamoreaux and his Los 
Alamos group jlU] developed a solid state technique that can improve the sensitivity to the electron 
EDM by 10^ to reach 10^^^ e cm. By operating at a lower temperature it is feasible to eventually 
reach a sensitivity of 10~^^ e cm, an improvement of eight orders of magnitude over the present 
sensitivity. The time scale for these is uncertain, as it is tied to funding prospects. Semertzidis and 
his Brookhaven group |U plan to trap muons in storage rings and increase the sensitivity of their 
EDM measurement by five orders of magnitude. A new measurement has been presented by the 
Sussex group [121 ■ ^ number of other experiments aim for an improvement in sensitivity by one or 
two orders of magnitude, and involve nuclear EDMs. 
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In order to understand the current and future constraints on the model, we can study what is 
the induced electron EDM, once we have satisfied the constraint from baryogenesis in ea. (|8c{|) : 

gg'{g^-g'^)kk'i^^lQ-^ (93) 
^ c 

The way we proceed is as follows. First, we fix g' = This is a good representative of the 

possible ratios between g and g\ as it is quite far from the region where the approximate symmetry 
in the case g = g' suppresses baryogenesis, and g' is not too small to suppress baryogenesis on its 
own. Later on we shall relax this condition. Having done this, we invert eq. (|93|) . to get an expression 
for g in terms of the other parameters of the model. Now, the couplings g, g' are expressed in terms 
of Tc, and in terms of kS and k' S which, as it will be useful, in the case of small mixing, can be 
thought of as respectively the singlet and the doublet mass. We shall impose the constraint on the 
couplings g,g',k,k' to be less than ~ 1, in order for these Yukawa couplings not to hit a Landau 
pole before the unification scale ~ 10^^ GeV. 

We decide to restrict our analysis to the case in which the possible ratios between the marginal 
couplings of the same sector of the theory are smaller than 2 orders of magnitude. Here, by same 
sector of the theory we mean either the CP violating sector, or the sector of the scalar potential 
and the strong gauge group. The justification of this relies in the fact that we wish to explore 
the most natural region of the parameter space. For this reason, we expect that there is no large 
hierarchy between the marginal couplings of the same sector. We think that 2 order of magnitudes 
is a threshold large enough to delimitate this natural region. However, since marginal couplings 
are radiatively relatively stable, and large hierarchies among them does not give rise to fine tuning 
issues, in principle, large hierarchies among the marginal couplings are acceptable. We think that, 
however, such a hierarchy would require the addition of further structure to the model to justify its 
presence, and we decide to restrict to the simplest realization of the model. As a consequence of 
this, the most important restrictions we apply are: 10~^ < k/k' < 10^, and 10^^ ^ f^s/^ ^ 10^- In 
particular, using ea. (|55|) and ea. H34() . this implies ^ Tc < S < 5 T^, which, using the limit k, k' < 1, 
implies that the largest of kS and k' S must be smaller than 5 Tc- 

The upper bound on the CP violating sector particles tells us that there is a small part of the 
parameter space which we are going to explore, in which the particles responsible for the production 
of baryons are non relativistic. In that case, we can approximately extend the result found in eq. (|93j) . 
with the purpose of having order of magnitude estimates, in the following way. From ea. H50() . it is 
clear that baryon production is Boltzmann suppressed if the particles are non relativistic. In that 
case, the CP violating charge will be generated by the scattering of these particle in the region where 
the induced mass on the particles is of the order of the critical temperature, as this is the condition 
of maximum CP violating interaction compatible with not being Boltzmann suppressed. Having 
observed this, it is easy to approximately extend the result of eq. (|93jl to the non relativistic case, by 
taking the vevs of the fields at the a value such that the induced mass is of the order of the critical 
temperature. 

In fig. 3, we show the induced EDM as a function of kS, for several values of the critical tem- 
perature Tc, for k' S = 100, on top, and k' S = 500, at the bottom, with the couplings g,g' chosen 
as explained above, in order to fulfil the baryogenesis requirement. Notice that 500 GeV is roughly 
the limit that LHC will put on SU(2) doublets. We choose the maximum CP violating phase. The 
horizontal lines represent the present constraint on EDM d(. < 1.7 x 10~'^'^e cm at 95% CL |37) . and 
the future expected one 1^1111011^ of order 10~^^e cm. A few features are worth to be noted. We 
see that, for fixed Tc, the EDM decreases as we decrease kS, for kS light enough. This is due the 
fact that, reducing kS, we reduce both the EDM and the produced quantity of baryons. However, 
the loss in the production of baryons is compensated with a much smaller, compared to the decrease 
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in kS, increase in the couplings g,g', so that, the baryon abundance can remain constant, while the 
EDM decreases. For large kS, the EDM decreases both because the mass of the particles in the 
loops becomes heavier and heavier, and also because the mixing becomes more and more suppressed. 
The maximum is located at the point where kS ~ k'S. In that case, in fact, the mass matrix has a 
diagonal piece roughly proportional to the identity, so, even though the diagonal elements are much 
larger than the off diagonal ones, mixing is much enhanced, and so is the EDM, which is proportional 
to the mixing. The lower and upper limit on the value of kS, as well as the minimum temperature, 
are dictated by the restrictions k'S < 5 Tc, fc5 < 5 Tc, and k/k' > 10~^. Now, as we decrease the 
critical temperature, the resulting EDM tends to decrease. In fact, decreasing the temperature much 
enhances the baryon production. This allows to decrease the couplings g,g' , which explains why the 
induced EDM decreases. We verify that increasing the hierarchy between g and g' does not change 
the results a lot, as the decrease in baryons production requires to make the other couplings large. 
Increasing the value of k'S up to the maximum allowed value of 1.1 TeV (since T^^^ ~ v) does not 
produce any relevant change in the result, because this raises the minimum temperature, so that 
baryogenesis requires larger couplings, which forces the EDM not to decrease relevantly. It is only 
once the CP violating phase is lowered to less than 10~^ that a small experimentally unreachable 
region is opened up around k'S = 500 GeV and kS < 1 GeV. The same region is opened also enlarg- 
ing the allowed hierarchy between the couplings to above 5 x 10^. However, clearly, the presence at 
the same time of a large hierarchy and of a very small CP phase requires some additional structure 
on the model to explain the reason for their presence. 

The main conclusion we can draw from combining the analysis on baryogenesis and EDM is 
that, at present, the most natural region of the parameter space is perfectly allowed, however, 
improvements in the determination of the EDM, are going to explore the entire viable region of the 
parameter space, so that absence of signal, would result in effectively ruling out the model, at least 
if no further structure is added. 

5 Comments on Dark Matter and Gauge Coupling Uni- 
fication 

The proposed model finds its main motivation in stabilizing the weak scale through the requirement 
of attaining baryons in our universe. However, further than this, it is clear that the model provides 
two other interesting phenomenological aspects: gauge coupling unification and dark matter. Here, 
we just briefly introduce these aspects, and we postpone a more detailed discussion to future work. 

Thanks to the two doublets we have inserted in our model, which have the same quantum numbers 
as higgsinos in the MSSM, gauge coupling unification works much better than in the standard model. 
As it was shown in gauge coupling unification with the standard model plus higgsinos works 
roughly as well as the MSSM at two loops level, with the only possible problem being the fact that 
the unification scale is a bit low at around ~ 10^^ GeV. The problem from proton decay can however 
be avoided with some particular model at the GUT scale US] . We expect that 2-loop gauge 
coupling unification works quite well also in this model, with only small corrections coming from the 
presence of the singlet scalar S. 

Concerning the Dark Matter relic abundance, the lightest of the newly introduced particles is 
stable, and so it provides a natural candidate for Dark Matter. Estimating the relic abundance 
is quite complex, as it depends on the composition of the particle, and on its annihilation and 
coannihilations rate. However, if our newly introduced Yukawa couplings are not very close to their 
upper limit (of order one), and if the lightest particle is mostly composed of the two doublets, then 
its relic abundance is very similar to the one of pure higgsino dark matter in split supersymmetry 
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Figure 3: Top, the predicted EDM given the constraint from Baryogenesis satisfied, as a function 
of kS for k' S = 100 GeV, g' = g/\/2, maximum CP violating phase, and Tc = 200 GeV (solid 
line), 100 GeV (long dashed), 25 GeV (short dashed). Bottom, the same for k' S = 500 GeV, and 
Tc = 200 GeV (sohd hne), 150 GeV (long dashed), 110 GeV (short dashed). The horizontal lines 
represent the present electron EDM constraint de < 1.7x 10~^^e cm at 95% CL |37j, and the expected 
improvement up to df. < 10^^ e cm [^[lUll^ • 



O Uni ISl IS im lini ) requiring the doublets to be around the TeV scale. In the case of singlet dark 
matter, we expect the singlet to give the right abundance for a much lighter mass, as it is naturally 
much less interacting. However, estimates become much more difficult, and we postpone the precise 
determination of the relic abundance to future work. 



6 Conclusions 

In this paper, we have addressed the solution to the electroweak hierarchy problem in the context 
of the landscape, following a recent model proposed in [TT] . 

We have shown that it is possible to connect the electroweak scale to a hierarchically small scale 
at which a gauge group becomes strong by dimensional transmutation. The assumption is that 
we are in a "friendly neighborhood" of the landscape in which only the relevant parameters of the 
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low energy theory are effectively scanned. We then realize the model in such a way that there is a 
fragile, though necessary, feature of the universe which needs to be realized in our universe in order 
to sustain any sort of life. As a natural continuation of Weinberg's "structure principle", the fragile 
feature is the presence of baryons in the universe, which are a necessary ingredient for the formation 
of clumped structures. We assume that, in the friendly neighborhood of the landscape in which we 
should be, baryogenesis is possible only through the mechanism of electroweak baryogenesis. Then, 
in order to produce a first order phase transition strong enough so that sphalerons do not wash 
out the produced baryon density, we develop a new mechanism to implement the electroweak phase 
transition. We introduce a new gauge sector which becomes strong at an exponentially small scale 
through dimensional transmutation. We couple this new sector to a singlet S which is then coupled 
to the higgs field. The electroweak phase transition occurs as the new gauge sector becomes strong, 
and produces a chiral condensation. This triggers a phase transition for the singlet S, which then 
triggers the phase transition for the higgs fields. In order to preserve baryon number, we need the 
phase transition to be strong enough, and this is true only if the higgs mass is comparable to the 
QCD scale of the strong sector. This solves the hierarchy problem. In order to provide the necessary 
CP violation, we introduce 2 SU(2) doublets ^-t with hypercharge ±1/2, and a gauge singlet s. 

When we require the model to describe our world, the model leads to falsifiable predictions. 

The main result of the paper is the computation of the produced baryon number, for which we 
obtain: 




The requirement that this baryon abundance should cope with the observed one leads to a lower 
bound on a combination of the product of the CP violating phase, the new couplings, the S vev, 
and the critical temperature Tc- 

We infer that Gauge Coupling Unification, and the right amount of Dark Matter relic abundance 
are easily achieved in this model. 

We study in detail the induced electron Electric Dipole Moment (EDM), and we find that, at 
present, the most natural region of the parameter space of the theory is allowed. However, soon 
in the future improvement in the EDM experiments will be sensitive the entire viable region of the 
parameter space of the model, so that absence of a signal would result in practically ruling out the 
model. 
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A Baryogenesis in the Large Velocity Approximation 

The method we have used in the main part of this work is based on some approximations that 
fail in limit of very fast wall speed: ~ 1. Even though this regime seems to be disfavored by 
actual computations of the wall speed|3T]. it is worth to try to estimate the result even in this 
case. In this appendix we are going to do an approximate computation in the large wall velocity 
approximation. In this new regime, calculations become much more complicated, as the assumption 
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of local thermodynamical equilibrium begins to fail, and the baryon number tends to be produced 
in the region of the wall, in the so called local baryogenesis scenario. 

We follow the treatment of In this approach, the phase transition is treated from an effective 
field theory point of view. A good way of looking at the high temperature phase of the unbroken 
phase is imagining it as discretized in a lot of cells, whose side is given by the typical size of the 
weak sphaleron barrier crossing configuration of the gauge field: ( ~ {awTc)~^ ■ Concentrating 
on each cell, we have that the thermal energy in a cell is in general much larger than the energy 
necessary to create a gauge field oscillation capable of crossing the barrier, and this means that most 
of the energy is in oscillations with smaller wavelength than So, these configurations cross the 
barrier at energy far above the one of the sphaleron configuration, and so their rate has nothing 
to do with the sphaleron rate, and this explain why their rate per unit volume is of order (In 
the case m| < 0, electroweak symmetry breaking has already occurred outside the bubble, but this 
discussion still roughly applies). We can parameterize the configuration in one cell with one variable 
T, which depends only on time. Obviously, this is a very rough approximation, but this is at least 
a beginning. We try to describe the dynamics of the configuration near crossing the barrier, at a 
maximum of the energy, that we fixe to be at r = 0. Near this point, we can write the following 
Lagrangian, which is very similar to the one in j51j : 



L{r,f) = ^f' + ^r' + ^bcPf (95) 



where Cj are dimensionless parameters depending on the different possible barrier crossing trajec- 
tories, while b will be determined shortly. Since the point r = is a maximum of the energy, no 
odd powers of r can appear. The appearance of the odd power in f can be understood because, at 
one loop level, the CP violating mass matrix of the particles in the CP violating sector induces an 
operator which contains the term TiF^uF'^'^ . This can be seen in the following way, in an argument 
similar to the one shown in |17] . We can imagine to do a chiral rotation on these field of amplitude 
equal to the CP violating phase. Because of the anomaly, this will induce the operator: 

O = ^(^TVF^.F^'^ (96) 
where 52 is the SU(2) weak coupling, and F^i, is the SU(2) field strength, and where 

= Arg(Det(M)) = Arcsm ( ,,.g2,,,(,). + + fefe-S^cos(^))^ ) ^''^ 

Promoting the vevs of h and S to the actual fields, we find the operator we were looking for. Now, 
the term TiF^yF^^ contains a term proportional to the time derivative of the Chern Simons number, 
and so the operator O must contribute to the effective action (|95|) with a term proportional to b(t)f, 

with b = Yg^, explaining the reason for the presence of the term in f 
The equation of motion for r is: 

'^--T7^'t> (98) 



If the wall is fast enough, we can solve it in the impulse approximation, to get: 

= --1^2^^ (99) 

This kick to At makes the distribution of velocities of barrier crossing configurations asymmetric, 
leading to a production of baryons with respect to antibaryons. 
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Now, this kick will be very inefficient in changing the distribution of baryons, unless the kick 
is larger than the typical speed tq a generic configuration would have if it crossed the barrier in 
the absence of the wall. So, we require: Af > tq. The fraction of configurations which satisfy this 
requirement is proportional to Af, but is very difficult to estimate. Following we just say that 
it is equal to /Af, where / is our "ignorance" coefficient. So, we finally get: 

UB ~ /AfC"^ (100) 

where we have reabsorbed the constants Cj into /. We finally get: 

There are a lot of heavy approximations which suggest that this estimate is very rough |51j : from 
the value of the coefficient Cj, to the approximation of restricting to one degree of freedom, to the 
impulse approximation in solving the differential equation, and to the estimate of the fraction of 
configurations influenced by the kick. All these suggests that we should take eq. (|l()lj) with / ~ 1, 
as at most an upper limit on the baryon production, as the authors of [SJ suggest. 
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